Whole-genome microarrays of fission yeast: characteristics, accuracy, reproducibility, and processing of array data

Background: The genome of the fission yeast Schizosaccharomyces pombe has recently been sequenced, setting the stage for the post-genomic era of this increasingly popular model organism. We have built fission yeast microarrays, optimised protocols to improve array performance, and carried out experiments to assess various characteristics of microarrays.|Results: We designed PCR primers to amplify specific probes (180-500 bp) for all known and predicted fission yeast genes, which are printed in duplicate onto separate regions of glass slides together with control elements (similar to13,000 spots/slide). Fluorescence signal intensities depended on the size and intragenic position of the array elements, whereas the signal ratios were largely independent of element properties. Only the coding strand is covalently linked to the slides, and our array elements can discriminate transcriptional direction. The microarrays can distinguish sequences with up to 70% identity, above which cross-hybridisation contributes to the signal intensity. We tested the accuracy of signal ratios and measured the reproducibility of array data caused by biological and technical factors. Because the technical variability is lower, it is best to use samples prepared from independent biological experiments to obtain repeated measurements with swapping of fluorochromes to prevent dye bias. We also developed a script that discards unreliable data and performs a normalization to correct spatial artefacts.|Conclusions: This paper provides data for several microarray properties that are rarely measured. The results define critical parameters for microarray design and experiments and provide a framework to optimise and interpret array data. Our arrays give reproducible and accurate expression ratios with high sensitivity. The scripts for primer design and initial data processing as well as primer sequences and detailed protocols are available from our website

Large-Scale Whole-Genome Sequencing Reveals the Genetic Architecture of Primary Membranoproliferative GN and C3 Glomerulopathy

BACKGROUND: Primary membranoproliferative GN, including complement 3 (C3) glomerulopathy, is a rare, untreatable kidney disease characterized by glomerular complement deposition. Complement gene mutations can cause familial C3 glomerulopathy, and studies have reported rare variants in complement genes in nonfamilial primary membranoproliferative GN. METHODS: We analyzed whole-genome sequence data from 165 primary membranoproliferative GN cases and 10,250 individuals without the condition (controls) as part of the National Institutes of Health Research BioResource–Rare Diseases Study. We examined copy number, rare, and common variants. RESULTS: Our analysis included 146 primary membranoproliferative GN cases and 6442 controls who were unrelated and of European ancestry. We observed no significant enrichment of rare variants in candidate genes (genes encoding components of the complement alternative pathway and other genes associated with the related disease atypical hemolytic uremic syndrome; 6.8% in cases versus 5.9% in controls) or exome-wide. However, a significant common variant locus was identified at 6p21.32 (rs35406322) (P=3.29×10−8; odds ratio [OR], 1.93; 95% confidence interval [95% CI], 1.53 to 2.44), overlapping the HLA locus. Imputation of HLA types mapped this signal to a haplotype incorporating DQA1*05:01, DQB1*02:01, and DRB1*03:01 (P=1.21×10−8; OR, 2.19; 95% CI, 1.66 to 2.89). This finding was replicated by analysis of HLA serotypes in 338 individuals with membranoproliferative GN and 15,614 individuals with nonimmune renal failure. CONCLUSIONS: We found that HLA type, but not rare complement gene variation, is associated with primary membranoproliferative GN. These findings challenge the paradigm of complement gene mutations typically causing primary membranoproliferative GN and implicate an underlying autoimmune mechanism in most cases

NMR Derived Model of GTPase Effector Domain (GED) Self Association: Relevance to Dynamin Assembly

Self-association of dynamin to form spiral structures around lipidic vesicles during endocytosis is largely mediated by its ‘coiled coil’ GTPase Effector Domain (GED), which, in vitro, self-associates into huge helical assemblies. Residue-level structural characterizations of these assemblies and understanding the process of association have remained a challenge. It is also impossible to get folded monomers in the solution phase. In this context, we have developed here a strategy to probe the self-association of GED by first dissociating the assembly using Dimethyl Sulfoxide (DMSO) and then systematically monitoring the refolding into helix and concomitant re-association using NMR spectroscopy, as DMSO concentration is progressively reduced. The short segment, Arg109 - Met116, acts as the nucleation site for helix formation and self-association. Hydrophobic and complementary charge interactions on the surfaces drive self-association, as the helices elongate in both the directions resulting in an antiparallel stack. A small N-terminal segment remains floppy in the assembly. Following these and other published results on inter-domain interactions, we have proposed a plausible mode of dynamin self assembly

Predicting Cellular Growth from Gene Expression Signatures

Maintaining balanced growth in a changing environment is a fundamental systems-level challenge for cellular physiology, particularly in microorganisms. While the complete set of regulatory and functional pathways supporting growth and cellular proliferation are not yet known, portions of them are well understood. In particular, cellular proliferation is governed by mechanisms that are highly conserved from unicellular to multicellular organisms, and the disruption of these processes in metazoans is a major factor in the development of cancer. In this paper, we develop statistical methodology to identify quantitative aspects of the regulatory mechanisms underlying cellular proliferation in Saccharomyces cerevisiae. We find that the expression levels of a small set of genes can be exploited to predict the instantaneous growth rate of any cellular culture with high accuracy. The predictions obtained in this fashion are robust to changing biological conditions, experimental methods, and technological platforms. The proposed model is also effective in predicting growth rates for the related yeast Saccharomyces bayanus and the highly diverged yeast Schizosaccharomyces pombe, suggesting that the underlying regulatory signature is conserved across a wide range of unicellular evolution. We investigate the biological significance of the gene expression signature that the predictions are based upon from multiple perspectives: by perturbing the regulatory network through the Ras/PKA pathway, observing strong upregulation of growth rate even in the absence of appropriate nutrients, and discovering putative transcription factor binding sites, observing enrichment in growth-correlated genes. More broadly, the proposed methodology enables biological insights about growth at an instantaneous time scale, inaccessible by direct experimental methods. Data and tools enabling others to apply our methods are available at http://function.princeton.edu/growthrate

The Princeton Protein Orthology Database (P-POD): A Comparative Genomics Analysis Tool for Biologists

Many biological databases that provide comparative genomics information and tools are now available on the internet. While certainly quite useful, to our knowledge none of the existing databases combine results from multiple comparative genomics methods with manually curated information from the literature. Here we describe the Princeton Protein Orthology Database (P-POD, http://ortholog.princeton.edu), a user-friendly database system that allows users to find and visualize the phylogenetic relationships among predicted orthologs (based on the OrthoMCL method) to a query gene from any of eight eukaryotic organisms, and to see the orthologs in a wider evolutionary context (based on the Jaccard clustering method). In addition to the phylogenetic information, the database contains experimental results manually collected from the literature that can be compared to the computational analyses, as well as links to relevant human disease and gene information via the OMIM, model organism, and sequence databases. Our aim is for the P-POD resource to be extremely useful to typical experimental biologists wanting to learn more about the evolutionary context of their favorite genes. P-POD is based on the commonly used Generic Model Organism Database (GMOD) schema and can be downloaded in its entirety for installation on one's own system. Thus, bioinformaticians and software developers may also find P-POD useful because they can use the P-POD database infrastructure when developing their own comparative genomics resources and database tools

De Novo Truncating Mutations in WASF1 Cause Intellectual Disability with Seizures.

Next-generation sequencing has been invaluable in the elucidation of the genetic etiology of many subtypes of intellectual disability in recent years. Here, using exome sequencing and whole-genome sequencing, we identified three de novo truncating mutations in WAS protein family member 1 (WASF1) in five unrelated individuals with moderate to profound intellectual disability with autistic features and seizures. WASF1, also known as WAVE1, is part of the WAVE complex and acts as a mediator between Rac-GTPase and actin to induce actin polymerization. The three mutations connected by Matchmaker Exchange were c.1516C>T (p.Arg506Ter), which occurs in three unrelated individuals, c.1558C>T (p.Gln520Ter), and c.1482delinsGCCAGG (p.Ile494MetfsTer23). All three variants are predicted to partially or fully disrupt the C-terminal actin-binding WCA domain. Functional studies using fibroblast cells from two affected individuals with the c.1516C>T mutation showed a truncated WASF1 and a defect in actin remodeling. This study provides evidence that de novo heterozygous mutations in WASF1 cause a rare form of intellectual disability

A gain-of-function variant in DIAPH1 causes dominant macrothrombocytopenia and hearing loss

Macrothrombocytopenia (MTP) is a heterogeneous group of disorders characterized by enlarged and reduced numbers of circulating platelets, sometimes resulting in abnormal bleeding. In most MTP, this phenotype arises because of altered regulation of platelet formation from megakaryocytes (MKs). We report the identification of DIAPH1, which encodes the Rho-effector diaphanous-related formin 1 (DIAPH1), as a candidate gene for MTP using exome sequencing, ontological phenotyping, and similarity regression. We describe 2 unrelated pedigrees with MTP and sensorineural hearing loss that segregate with a DIAPH1 R1213* variant predicting partial truncation of the DIAPH1 diaphanous autoregulatory domain. The R1213* variant was linked to reduced proplatelet formation from cultured MKs, cell clustering, and abnormal cortical filamentous actin. Similarly, in platelets, there was increased filamentous actin and stable microtubules, indicating constitutive activation of DIAPH1. Overexpression of DIAPH1 R1213* in cells reproduced the cytoskeletal alterations found in platelets. Our description of a novel disorder of platelet formation and hearing loss extends the repertoire of DIAPH1-related disease and provides new insight into the autoregulation of DIAPH1 activity.The NIHR BioResource- Rare Diseases and the associated BRIDGE genome sequencing projects are supported by the National Institute for Health Research (NIHR; http://www.nihr.ac.uk). B.N. was supported by the Deutsche Forschungsgemeinschaft (SFB 688). S.S. was supported by a grant of the German Excellence Initiative to the Graduate School of Life Sciences, University of Würzburg. ET, DG, JCS, SP, IS, CJP, RM, SAsh, ST and KS are supported by the NIHR BioResource - Rare Diseases. KF, CT, and CVG are supported by the Fund for Scientific Research-Flanders (FWO-Vlaanderen, Belgium, G.0B17.13N) and by the Research Council of the University of Leuven (BOF KU Leuven‚ Belgium, OT/14/098). WNE is supported by the Cancer Council Western Australia. Research in the Ouwehand laboratory is supported by program grants from the European Commission, NIHR to WJA, SM, MK, RP, SBJ and WHO under numbers RP-PG-0310-1002; the laboratory also receives funding from NHS Blood and Transplant; CL and SKW are supported by Medical Research Council (MRC) Clinical Training Fellowships (number MR/K023489/1) and TKB by a British Society of Haematology/NHS Blood and Transplant grant. MAL and CL are supported by the Imperial College London Biomedical Research Centre; JRB acknowledges support by the NIHR Cambridge Biomedical Research Centre and SR by the Medical Research Council and Cambridge Biomedical Research Centre. CVG is holder of the Bayer and Norbert Heimburger (CSL Behring) Chairs. ADM is supported by the NIHR Bristol Cardiovascular Biomedical Research Unit

Platelet function is modified by common sequence variation in megakaryocyte super enhancers

Linking non-coding genetic variants associated with the risk of diseases or disease-relevant traits to target genes is a crucial step to realize GWAS potential in the introduction of precision medicine. Here we set out to determine the mechanisms underpinning variant association with platelet quantitative traits using cell type-matched epigenomic data and promoter long-range interactions. We identify potential regulatory functions for 423 of 565 (75%) non-coding variants associated with platelet traits and we demonstrate, through ex vivo and proof of principle genome editing validation, that variants in super enhancers play an important role in controlling archetypical platelet functions

Whole genome sequencing reveals that genetic conditions are frequent in intensively ill children

Purpose With growing evidence that rare single gene disorders present in the neonatal period, there is a need for rapid, systematic, and comprehensive genomic diagnoses in ICUs to assist acute and long-term clinical decisions. This study aimed to identify genetic conditions in neonatal (NICU) and paediatric (PICU) intensive care populations. Methods We performed trio whole genome sequence (WGS) analysis on a prospective cohort of families recruited in NICU and PICU at a single site in the UK. We developed a research pipeline in collaboration with the National Health Service to deliver validated pertinent pathogenic findings within 2andndash;3andnbsp;weeks of recruitment. Results A total of 195 families had whole genome analysis performed (567 samples) and 21% received a molecular diagnosis for the underlying genetic condition in the child. The phenotypic description of the child was a poor predictor of the gene identified in 90% of cases, arguing for gene agnostic testing in NICU/PICU. The diagnosis affected clinical management in more than 65% of cases (83% in neonates) including modification of treatments and care pathways and/or informing palliative care decisions. A 2andndash;3andnbsp;week turnaround was sufficient to impact most clinical decision-making. Conclusions The use of WGS in intensively ill children is acceptable and trio analysis facilitates diagnoses. A gene agnostic approach was effective in identifying an underlying genetic condition, with phenotypes and symptomatology being primarily used for data interpretation rather than gene selection. WGS analysis has the potential to be a first-line diagnostic tool for a subset of intensively ill children.</p

De Novo Truncating Mutations in WASF1 Cause Intellectual Disability with Seizures

© 2018 Next-generation sequencing has been invaluable in the elucidation of the genetic etiology of many subtypes of intellectual disability in recent years. Here, using exome sequencing and whole-genome sequencing, we identified three de novo truncating mutations in WAS protein family member 1 (WASF1) in five unrelated individuals with moderate to profound intellectual disability with autistic features and seizures. WASF1, also known as WAVE1, is part of the WAVE complex and acts as a mediator between Rac-GTPase and actin to induce actin polymerization. The three mutations connected by Matchmaker Exchange were c.1516C>T (p.Arg506Ter), which occurs in three unrelated individuals, c.1558C>T (p.Gln520Ter), and c.1482delinsGCCAGG (p.Ile494MetfsTer23). All three variants are predicted to partially or fully disrupt the C-terminal actin-binding WCA domain. Functional studies using fibroblast cells from two affected individuals with the c.1516C>T mutation showed a truncated WASF1 and a defect in actin remodeling. This study provides evidence that de novo heterozygous mutations in WASF1 cause a rare form of intellectual disability